Automatic Control and System Engineering Journal, ISSN 1687-4811, Volume 17, Issue 1, ICGST, Delaware, USA, June 2017 



ica«T 


www.icgst.com 


ACi? = 


Digital Optimal Compensation of Process Control Systems 


Kamen M. Yanev 1 

Department of Electrical Engineering, University of Botswana, Gaborone, Botswana 
yanevkm@yahoo.com, yanevkm@mopipi.ub.bw 
http://www.ub.bw 


Abstract 


1. Introduction 


This paper is suggesting a practical method for realizing 
optimal digital compensation of process control systems. 
Initially, a cascade continuous compensator prototype is 
designed and further it is converted into its digital 
equivalent. The purpose of the compensator is to 
eliminate some properly selected poles of the original 
control system transfer function, to introduce new 
dominant poles and specific system amplification. As a 
result, the system's transient response and stability are 
considerably improved. The control system’s relative 
damping ratio, percent maximum overshoot, phase 
margin, rise and settling times, are optimized. 

Keywords: Marginal system; Dominant poles; Stability; 
Analogue prototype; Lag-Lead multistage compensation; 
Digital compensation; Optimal system; Transient 
responses; 

Nomenclature 

ITAE Integral of the Time- weighted Absolute Error 

£ (DR) Closed-loop relative damping ratio 

PMO Percent Maximum Overshoot 

t s Settling time 

t m Maximum overshoot time 

e ss Steady-State error 

Gpfs) Transfer function of the plant open-loop system 
s Laplace operator 

K, K' Variable gains of the linear system 
Gp(s) Characteristic equation of the open-loop system 
G c (s) Transfer function of the analogue compensator 
G c '(s) Transfer function of the multi-stage lead section 
of the compensator continuous prototype 
G c "(s) Transfer function of the attenuation 
G c ’"(s) Transfer function of the of the lag compensation 
and amplification 
T s Sampling period 

T min Minimum time-constant of the linear section of 

the system 

cos Sampling frequency 
z z-transform operator 

Gcd(z) Transfer function of the series compensation 
controller stage in the discrete-time domain 


High quality performance of control systems can be 
achieved by applying the suggested cascade optimal 
digital compensation. It is implemented by employing a 
number of steps for the design of a proper digital 
compensator. The purpose of the compensator is to 
eliminate some properly selected poles of the system’s 
transfer function. At the same time it introduces new 
dominant poles and a specifically designed amplification. 
As a result, the quality of the system's performance in 
terms of its transient response, stability and accuracy are 
considerably improved. The suggested technique of 
successive steps brings the system to its marginal state, 
further subjecting it to multi-stage phase-lead, phase-lag 
compensation and amplification. The continuous cascade 
compensator prototype is finally converted into its digital 
equivalent and implemented as a microcontroller 
compensator. The main purpose of the suggested method 
is the compensated control system to meet the ITAE 
criterion and accordingly, the following objectives [1]: 

C = 0.707; (PMO) < 4%; *,/*„< 1.4; e ss < 1 % (type 0) 

Phase-lead and phase-lag compensation technique is a 
well known mechanism for improving control systems’ 
performance as demonstrated in considerable number of 
sources [1], [2], [3], [4]. The disadvantage of this method 
of compensation, as described in the sources, is that it 
only recommends the range of the compensator’s 
frequencies and it depends on the designer’s expertise to 
achieve the best performance of the compensated system. 
Further, the direct application of the compensation 
technique causes compromise with the system’s gain. 

Another method to meet the objectives of the ITAE 
criterion is to apply a PI or a PID controller to process 
control systems. This is a well known approach and it is 
also described in the sources mentioned above. Different 
methodologies are recommended for determination of the 
PI or PID controller constants where results are usually 
achieved with a lot of assumptions. Finally, the best 
system performance is proved to be accomplished by a 
procedure of manual MATLAB interaction tuning of the 
PI or PID controller constants [5]. 
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The importance of the new controller design technique, 
presented by the author of this research, is that it is 
suggesting a straight forward and exact achievement of 
the ITAE criterion objectives. Further contribution of the 
suggested by the author technique is the conversion of 
the achieved controller from its continuous to its digital 
equivalent. In this way the compensation is implemented 
with the aid of microcontrollers that are programmed 
based of the difference equations, representing the digital 
transfer functions of the controller stages. The application 
of digital control is in accordance to the Euler's 
approximation [6], [7], [8]. It recommends that the 
sampling period T s of the discrete system, should be 
within the range T s < (0. 1 T min to 0.2 T min ), where T min is the 
minimum time-constant of the continuous-time system or 
the analogue plant model prototype. Only then there will 
be a close match between the discrete and the 
continuous-time system performance and both, the 
optimal plant analogue prototype and its optimal digital 
equivalent are having exactly the same relative damping 
of £= 0.707. If the Euler's approximation is applied, it is 
apparent that there is a complete match between the 
corresponding pole location on the s-plane and on the z- 
plane. Due to this conclusion, the controller stages are 
designed in the continuous-time domain and further 
converted into their discrete-time equivalents [9]. 

For the clarity of the presented research, this paper is 
organized in the following sequence. Initially, the 
compensation design for a system Type 0 is described for 
the cases of considerable difference and insignificant 
difference between dominant and insignificant system 
poles. With the aid of three successive steps the system’s 
performance is preliminary evaluated. Further, the 
implementation of five successive rules are bringing the 
system’s performance to perfection, satisfying the ITAE 
criterion objectives. This is followed by the application 
of the compensation technique to a system Type 1. 
Finally, a pseudo code technique is applied for describing 
the different system cases with the aid of a flowchart, in 
this way clarifying the paths and steps that are to be 
followed for the design of the compensation. 

2. Compensator Design for a System Type 0 

Case 1: Considerable Difference between Dominant and 
Insignificant Poles 

The compensation design can be demonstrated for a case 
of a cruise control system [10] with a transfer function: 

K ' 

Gp\ (s) = = 

(i + r^(i + rs)(i + rs) 

1 2 3 


(l + 0.D)(l + 0.2s)(l + s) 

500 _ K 

s 3 + 16s 2 + 65s + 50 s 3 + 16s 2 + 65s + 50 


Step 1: Evaluation of the Original Control System 

The system’s performance is unacceptable due to large 
oscillations as seen from Figure 1. The system’s transient 
response is achieved by applying the following 
MATLAB code: 



» Gpl=tf([0 500], [1 16 65 50]) 
» Gfbl=feedback(Gpl,l) 

» step(Gfbl) 


Step Response 



Time (sec) 

Figure 1: Transient Response of the Original Control System 

Step 2: Determination of the Marginal System Gain 

To apply the suggested compensation technique, the 
system initially has to be brought to its marginal state by 
tuning its original gain K\ The critical value of K' is 
achieved by means of applying the proposed by the 
author, in some of his previous publications, method of 
Advanced D-Partitioning [11], [12], [13]. 

To apply this method, the characteristic equation of the 
system (2) is derived from Equation (1), considering the 
system’s gain value K' = 0.02K as unknown. 

G P (s) = s 3 +I6s 2 +65s + 50 + K (2) 

From Equation (2) the unknown value of K can be 
presented as follows: 

K = -s 3 -16s 2 -65s -50 (3) 

The D-partitioning curve in terms of the variable 
parameter K can be plotted in the complex plane within 
the frequency range -oo < 00 < +oo facilitated by 
MATLAB the “nyquist” m-code, as seen in Figure 2. 


» K=tf([-1 -16 -65 -50], [0 1]) 
Transfer function: 

-s A 3 -16 s A 2 - 65 s - 50 
» nyquist(K) 


D-Partitioning in Terms of the Variable K 



Figure 2: D-Partitioning Facilitated by the “nyquist” m-code 
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To avoid any misinterpretation of the D-Partitioning 
procedure, the “nyquisf ’ m-code can be modified into a 
“dpartition” m-code with the aid of the MATLAB Editor 
and a proper formatting. 

It is seen that the D-partitioning determines three regions 
on the AT-plane: D(0), D(l) and D(2). Applying the 
Advanced D-Partitioning method, only D(0) is the region 
of stability, being always on the left-hand side of the 
curve for the frequency range from -oo to +oo [11], [12]. 

At K = 990, the system is marginal and the system’s gain 
becomes K' = 0.02 K = 0.02x990 = 19.8. 

Comparison between the original and the marginal 
system transient responses is shown in Figure 3. 

» Gpl=tf([0 500], [1 16 65 50]) 

» Gfbl=feedback(Gpl,l) 

» Gp2=tf([0 990], [1 16 65 50]) 

» Gfb2=feedback(Gp2,l) 

» step(Gfbl,Gfb2) 


Step Responses of the Original and Marginal Systems 



Figure 3: Comparison between the Original and the Marginal Linear 
Prototype Systems’ Transient Responses 


Considering Equation (1), the transfer function of the 
marginal control system is presented as: 


Gp2( s ) 


19.8 

(1 + 0.1s)(l + 0.2s)(l + s) 


( 4 ) 



3 1 ' 3 2 

DR ~ P MO ts/tm | 


Figure 4: Determination of the Optimum Values of ai and on 


The most dominant pole of the open-loop system should 
be left uncompensated. Since the plant transfer function 
from Equation (1) is of a third order, according to Rule 1, 
the multi-stage lead section of the compensator 
continuous prototype should have a transfer function: 

' (1 + a,T,s){\ + a ? T ? s) 

G (s) = - 2 -^~ = 

a a (1 + ZA)(1 + T^s) 

12 (5) 

(1 + O.lsXl + 0.2s) 

” 100(1 + 0.01^(1 + 0.02s) 


Rule 2: The current gain is maintained by amplification 
equal to the product aia20i3...[13], [14]. 

G\s) = aa =10x10 = 100 (6) 

c 1 2 


Rule 3: To optimize a marginal closed-loop system, 
single-stage section of the continuous prototype lag 
compensator with factor /? = 10 is to be applied for a 
zero-pole cancellation of the uncompensated most 
dominant pole of the open-loop system. 

To optimize the steady- state error e ss the current gain 
should be increased by the factor ^=10 (proved from 
Figure 5 [13], [14]. Both relationships shown in Figure 4 
and Figure 5 are achieved by tracking procedures. 


Step 3: Optimization of £ , t s /t m and PMO 


Initially, the cascade compensator is designed as a 
continuous prototype and further it is converted into its 
digital equivalent. The following rules are applied: 


Rule 1: To optimize £ , t s / t m and the PMO of a Type 0 
marginal closed-loop system, a cascade multi-stage 
compensation with factors of a 1 , 2 , 3 ,.. = 10 is applied for a 
zero-pole cancellation. 

This rule is considered from the graphical relationship of 
Figure 4 [13], [14]. 


The number of the compensating stages N is to be one 
less than the order of the open-loop system, i.e. N = n + 


m- 1 . 
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The most significant pole in Equation (1) is p 3 = -1. 
Applying Rule 2, the section of the lag compensation and 
amplification is presented by: 


H1 + 7V) - 10(1 + s') 
(i +pr/> (i + i Os) 


(7) 


The transfer function of the full compensator is the 
product of the transfer functions of the three sections: 


G (s) = G (s) x G (s) x G (s) = 

c c C c 

10(1 + O.Ls'Xl + 0.2s)(l + s) V y 

” (l + 0.0Ls)(l + 0.02s)(l + 10s) 

Applying the full compensation, taking into account 
equations (4) and (8), the transfer function of the open- 
loop system becomes: 


G(s) = G ( s)xG (s) = 

c p 

100 

(1 + 0 . 01 j )(1 + 0 . 02 j )(1 + 10 j ) 
100 


(9) 


0.002s 3 + 0.3002 s 2 + 10.03s + 1 

Comparison between the transient responses of the 
original, the marginal and the compensated system, as 
shown in Figure 6, is achieved by the code: 


» Gpl=tf([0 500], [1 16 65 50]) 
» Gfbl=feedback(Gpl,l) 
Transfer function: 

500 

s A 3 + 16 s A 2 + 65 s + 550 
» Gp2=tf([0 990], [1 16 65 50]) 
» Gfb2=feedback(Gp2,l) 
Transfer function: 

990 


s A 3 + 16 s A 2 + 65 s + 1040 
» Gp3=tf([0 100], [0.002 0.3002 10.03 1]) 
» Gfb3=feedback( Gp3,l) 

Transfer function: 

100 

0.002 s A 3 + 0.3002 s A 2 + 10.03 s + 101 
» step(Gfbl,Gfb2,Gfb3) 


Step Responses of the Original. Marginal and 
the Compensated System 



Figure 6: Comparison between the Transient Responses of the Original, 
Marginal and Compensated Linear Prototype Systems 


The Bode diagram of the compensated system is plotted 
by applying the code: 

» Gp3=tf([0 100], [0.002 0.3002 10.03 1]) 

» margin (Gp3) 

Bode Diagram 


Gm = 22 dB (at 70.8 rad/sec) , Pm = 70.8 deg (at 1 1 .6 rad/sec) 



Figure 7: Bode Diagrams of the Compensated System Proving 
Phase Margin PM = 70.8 and Damping £ « 0.707 


The system performance after the compensation, seen 
from the transient response in Figure 6 and the Bode 
diagram in Figure 7, demonstrates the achievement of 
optimal performance [1], [2], [14]. 


Step 4: Conversion of the Continuous Prototype 
Compensated System into its Digital Equivalent. 


In order to apply digital compensation control, in 
accordance to the Euler's approximation [15] the 
sampling period T s of the discrete system, should be 
within the range T s < (0. 1 T min to 0.2 T min ), where T min is the 
minimum time-constant of the continuous-time system. 

By applying the Euler's approximation, a close match is 
expected between the discrete and the continuous-time 
system performance. For the case under discussion, the 
plant’s minimum time constant is T mi n = 0.1 sec and the 
sampling period is chosen as T s = 0.01 sec, satisfying the 
condition T s < 0.1 T min . When analysis is applied in the 
discrete-time domain, the considered frequency variation 
is within the range co = ± coJ2 = ±2n/2T s [6], [16]. 


The analysis in the discrete-time domain is applied with 
the aid of the Bilinear Tustin Transform [16], [17] that is 
a first-order approximation of the natural logarithm 
function, being an exact mapping of the z-plane to the s- 
plane and vice versa. It maps positions from the jco axis 
in the s-plane to the unit circle |z| = 1 in the z-plane. The 
Bilinear Tustin approximation and its inverse are 
presented as [16]: 


S T s e 

z=e s = — 


\ + sTJ2 m 
l-sTJ2' 


An(z) = 

1 s 


2 _ 

T s 


z — 1 
z + 1 


( 10 ) 


To employ the Bilinear Tustin Transform, the MATLAB 
'tustin' code is to be applied, performing the transform of 
the continuous-time to discrete-time system equivalents. 
The results demonstrated in Figure 6 and Figure 7 can be 
compared with the results after converting the analogue 
system model prototype into its digital equivalent: 
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» Gpl=tf([0 500], [1 16 65 50]) 

» Gfbl=feedback(Gpl,l) 

» Gp3=tf([0 100], [0.002 0.3002 10.03 1]) 

» Gfbld = c2d(Gfbl,0.01, Austin ) 

Transfer function: 

5. 778e-005 z A 3 + 0. 0001 733 z A 2 + 0. 0001 733 z + 5. 778e-005 

z A 3 - 2.846 z A 2 + 2.698 z - 0.852 
Sampling time: 0.01 
» Gfb3d = c2d(Gfb3,0.01, Austin ) 

Transfer function: 

0.003321 z A 3 + 0.009962 z A 2 + 0.009962 z + 0.003321 

z A 3 - 1.916 z A 2 + 1.139 z - 0.1958 
Sampling time: 0.01 
» step(Gfbld,Gfb3d) 


Step Responses of the Original and the 
Compensated Digital System 



Figure 8: Comparison between the Transient Responses of the 
Original and Compensated Digital Systems 

The Bode diagram of the compensated digital system is 
plotted by applying the code: 

» Gp3=tf([0 100], [0.002 0.3002 10.03 1J) 

» Gp3d = c2d( Gp3, 0. 01,tustin ) 

Transfer function: 

0.003332 z A 3 + 0.009995 z A 2 + 0.009995 z + 0.003332 


z A 3 - 1.932 z A 2 + 1.132 z - 0.1998 
Sampling time: 0.01 
» margin (Gp3d) 


Bode Diagram 



Frequency (rad/sec) 


Figure 9: Bode Diagrams of the Digital Compensated System Proving 
Phase Margin PM = 70.8 and Damping C , « 0.707 


The results obtained from Figure 8 and Figure 9, 
compared with the results from Figure 6 and Figure 7, 
prove the close match between the performance of the 
continuous-time model prototype and the digital model of 
the system (PM = 70.8 and « 0.707). 

The performance specifications of the compensated 
system, determined from the transient response and the 
Bode diagram shown in Figure 8 and Figure 9, are 
compared with the objectives for optimal performance. 

Table 1 

Objectives and Real Results for Case 1 


Specifications 

Objectives 

Real Results 

Consideration 

c 

= 0.707 

= 0.707 

Matching 

PMO 

<4% 

= 2.8% 

Better 

ts(l%)/tm 

< 1.49 

= 1.25 

Better 

e ss (t) 

< 1% 

= 0.06% 

Better 


From the graphs at Figure 8 and Figure 9 and from the 
summary in Table 1, it is seen that the transient response 
of the compensated system in terms of relative damping 
ratio £ percent maximum overshoot (PMO), time ratio 
ts(i%)/t m and steady- state error e ss (t) is either matching or it 
is better than the one of the set objective. 

Step 5: Design of a Digital Compensation Controller 
based on a Microcontroller (Case 1). 

The design of a digital compensation controller founded 
on its analogue prototype is realized again by applying 
the Tustin bilinear transform [16], [17]. From the 
Equation (8), the transfer function of the modified 
compensation controller stage can be presented as: 

_ 10(l + 0.Ls)(l + 0.2,s)(l + ,s) 

C~ (l + 0.0Ls)(l + 0.02^(1 + 105) ” 

( 11 ) 

2s 3 + 32s 2 + 130.S + 100 

0.02s 3 + 3.002s 2 + 100.3s + 10 

The digital equivalent of the series compensator 
controller stage [16], [17] is realized by the code: 


» Gc=tf([2 32 130 100], [0.02 3.002 100.3 10]) 
» Gcd = c2d(Gc,0.1, f tustin ) 

Transfer function: 

9.328 z A 3 - 17.15 z A 2 + 9. 743 z - 1.688 

z A 3 + 0.1052 z A 2 - 0. 7986 z - 0.2829 
Sampling time: 0.1 


Following the code result, the transfer function of the 
series compensation controller stage in the discrete-time 
domain is as follows: 

G 9.328z 3 -17.15z 2 + 9.743z - 1.688 _ Y(z) ^ 2 ) 

z 3 +0.1052z 2 -0.7986z -0.2829 


To enable the implementation of a microcontroller that is 
connected in cascade with the original plant, the transfer 
function of the series digital robust control stage is 
represented by the following difference equation [18]: 
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y(kT) = 0.1052 y[{k - l)r]- 0.7986 y[{k - 2 )r] + 

- 0.2829 y[(k-w]+ (13) 

+ 9328x(kT) - \l.\5x[(k - l)r] + 

+ 9.743x[(/t - 2)T\- 1.688x[(/t - 3)r] 

Case 2: Insignificant Difference between Dominant and 
Insignificant Poles 

Rule 4: If the most significant pole of the open-loop 
system has a real part close in value to that of an 
insignificant pole, Rule 3 is modified to Rule 4, where 
successive two-stage lag compensation with factors pi = 
p 2 and a factor 8 should be applied. 

The values of the factors pi and p 2 at Rule 4 are chosen by 
the same considerations as in Rule 2. The value of the 
factor 8 is varied by a tracking procedure on the transient 
response, searching for the optimum performance of the 
compensated system. It is seen from Figure 10 that the set 
of objectives can be met if the factor <5*= 80 [13], [18]. 
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Figure 10: Determination of Optimum Value of 8 

Rule 4 can be illustrated for a system Type 0 [19] of an 
armature control dc motor and a type driving 
mechanism with original gain K = 10. Applying the D- 
partitioning technique it is brought to its marginal state, 
where its gain becomes K = 20, as demonstrated below: 

20 

G P3 (s) = (14) 

(1 + 0.01s)(l + 0.02s)(l + 0. Is) 

The real part of the most significant pole is p 3 = -10 and it 
is only 5 times smaller than the real part of the one of the 
insignificant poles pi = -50. 

Again, in accordance to Rule 1 , a two-stage lead 
compensation and attenuation is applied, where the factors 
are ai - (X2= 10. Following Rule 2, the gain is maintained 
by amplification equal to the product otiOt 2 . Further, 
following Rule 4, a two-stage lag compensation with 
factors pi = p 2 = 10 and a factor amplification £=80 is 
suggested. Using a related sequence and applying 
Equation (5), the two-stage lead stage is presented by: 


Applying Rule 2, the second compensation section is 
achieving the amplification: 

G(s) = aa =10x10 = 100 (16) 

c 1 2 

Applying Rule 4, the two-stage lag compensation and 
amplification is presented by: 

G '" () <y ( 1 + VX 1 + r /) 80(1 + 0.15X1 + *) (17) 

c W (1 + PT^s)(\ + PTs) (1 + *)(! + 10s) 

The transfer function of the complete compensator as a 
product of the three sections is: 

G (s) = G (s) x G(s) xg'\s) = 

80(1 + 0.01s)(l + 0.02s)(l + 0. Is) 

~ (l + 0.001s)(l + 0.002s)(l + 10s) 

After applying the full compensation, 

Equations (14) and (18), the analogue prototype of the 
transfer function of the open-loop compensated system 
becomes: 

G(s) = G (s) x G (s) = 

P 1600 (19) 

(1 + 0.00U)(1 + 0.002s)(l + 10s) 

In the same way the Bilinear Tustin Transform is applied, 
converting the continuous-time system prototype into its 
digital equivalent. The system’s transient response before 
and after compensation is shown in Figure 11 and is 
obtained by the following code: 

» Gpl=500000/(s+10)/(s+50)/(s+100) 

» Gpfbl=feedback(Gpl,l) 

» G=80000000/(s+0.1)/(s+500)/(s+1000) 

» Gfb=feedback(G,l) 

» Gpfbld = c2d(Gpfbl, 0.002, 'tustin ) 

» Gfbd = c2d(Gfb, 0.002, 'tustin ') 

» step(Gpfbld,Gfbd) 


(18) 

considering 


Step Responses of the Original and the 
Compensated Digital Control System 



' (1 + a,T,s)(l + a^T^s) 

GJs) = ^2- = 

a a (1 + 7^)(1 + T 2 s) 

(1 + 0 . 0 D )(1 + 0 . 02 s ) 


100(1 + 0.001s)(l + 0.002) 


Figure 1 1 : Comparison between the Transient Responses of the 
Original and Compensated Digital System (Case 2) 

The Bode diagram of the compensated digital system is 
shown at Figure 1 1 and plotted by applying the code: 
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»G=80000000/(s+0.1)/(s+500)/(s+1000) 
» Gd = c2d(G, 0.002, Austin ) 
Zero/pole/gain : 

0.026664 (z+l) A 3 


z (z-1) (z-0.3333) 
Sampling time : 0. 002 


Bode Diagram 


Gm = 1 9.4 dB (at 61 6 rad/sec) , Pm = 64.6 deg (at 1 50 rad/sec) 



Frequency (rad/sec) 


Figure 12: Bode Diagrams of the Digital Compensated System Proving 
Phase Margin PM = 64.6 and Damping £ « 0.646 

From the summary in Table 2, it is obvious that the 
objectives are met. Again, the real results for the 
compensated control system are either close or better 
than the set specifications. 

Table 2 

Objectives and Real Results for Case 2 


Specifications 

Objectives 

Real Results 

Consideration 

c 

= 0.707 

= 0.646 

Close 

PMO 

<4% 

= 2.8% 

Better 

ts(l%)/tm 

< 1.49 

= 1.25 

Better 

e ss (t) 

<1% 

= 0.06% 

Better 


Step 5: Design of a Digital Compensation Controller 
based on a Microcontroller (Case 2). 


The digital equivalent of the series compensator 
controller stage can be realized in a similar procedure, 
like the one implemented for example of Case 1. The 
transfer function compensation controller stage sown in 
equation (18) is modified and can be presented in the 
following format: 

_ 80(1 + 0.01s)(l + 0.02s)(l + 0.1s) _ 
c~ (1 + 0.00Ls)(l + 0.002.s)(l + I0s) ~ 

( 20 ) 

_ 0.0016s 3 + 0.256s 2 + 10.% + 80 

0.00002s 3 + 0.030002 s 2 + 10.003s + 1 


Considering the equation (20), the digital equivalent of 
the series compensator controller stage is realized by the 
code: 


» Gc=tf([0.0016 0.256 10.4 80], [0.00002 0.030002 10.003 1J) 
Transfer function: 

0.0016 s A 3 + 0.256 s A 2 + 10.4 s + 80 

2e-005 s A 3 + 0.03 s A 2 + 10s + l 
» Gcd = c2d(Gc,0.01, Austin ) 

Transfer function: 

1.891 z A 3 + 1.441 z A 2 - 0.1501 z - 0.1801 


z A 3 + 0.8938 z A 2 - 0.9782 z - 0.8781 
Sampling time: 0.01 

Following the code result, the transfer function of the 
series compensation controller stage in the discrete-time 
domain is as follows: 


1.891z 3 + 1.441z 2 - 0.1501 z - 0.1801 Y(z) n\) 

G cd (A = = f 1 

z 3 + 0.8938z 2 - 0.9782z - 0.8781 

Taking into accounting equation (21), the transfer 
function of the series digital robust control stage is 
represented by the following difference equation (22). A 
microcontroller, connected in cascade with the original 
plant, can be programmed in accordance with this 
difference equation [18], [20]: 

y(kT) = 0.8938 y[{k - l)j]- 0.9782 y[(k - 2 )r] + 

- 0.8781 >>[(& - 3)J ]+ ^2) 

+ 1.891x(*r) + 1.44l4(/fc-l)r] + 

- 0.1501x[(fc - 2)T]- 0.1801 x[(k - 3)j] 


3. Compensator Design for a System Type 1 


The compensation technique for a light tracking system 
with a transfer function of Type 1 can be demonstrated 
with the aid of the following transfer function [21], [22]: 
K' 50 

G pa (s) = = = 

5(1 + 7A)(1 + Ts) 5(1 + 0.025)(1 + 0.055) ^3) 

500 K 

S 3 +7s 2 +1005 s 3 +7s 2 +100s 

Step 1 : Evaluation of the Original Control System 

The demonstration of the system’s transient response is 
achieved by applying the following code: 


»Gp4=tf([0 500], [1 7100 0]) 
Transfer function: 

500 


s A 3 + 7 s A 2 + 100 s 
» Gfb4=feedback(Gp4,l) 
»step(Gfb4) 


Step Response 



Time (sec) 


Figure 13: Transient Response of the Original System Type 1 
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Step 2: Determination of the Marginal System Gain 

The critical marginal value of the unknown parameter K 
is achieved by means of applying the method of 
Advanced D-Partitioning. From equation (23), the 
characteristic equation of the system is obtained as: 

G / ,(5) = s 3 +7s 2 +1005 + /: (24) 

Reflecting on equation (24) the unknown value of K can 
be presented as follows: 

K = -s 3 -Is 2 -1005 (25) 

The D-partitioning curve in terms of the variable 
parameter K can be plotted in the complex plane within 
the frequency range -oo < co < +oo facilitated by the 
following code and seen in Figure 14. 


» K=tf([-1 -7 -100 0],[0 1]) 
Transfer function: 

-s A 3 - 7 s A 2 - 100 s 
» nyquist(K) 


D-Partitioning in Terms of the Variable K 



Figure 14: D-Partitioning in Terms of the Parameter K 

The D-partitioning determines three regions on the K- 
plane: D(0), D(l) and D(2). Only D(0) is the region of 
stability, being the one, always on the left-hand side of 
the curve for a frequency variation from -oo to +oo. 

At K = 700, the system is marginal and the system’s gain 
becomes K' = 0. IK = 0. 1 x700 = 70. 


Then equation (23) of the plant transfer function of Type 
1 is modified to its marginal case as follows: 


C?) 


r 


s(l + Ts)(l + T s) 
1 2 
70 


(26) 


Rule 5:_The pure integration or the most significant pole 
of the open-loop system should be left uncompensated. 
The existing system gain should be maintained by 
amplification equal to the product sa 1012013..., where the 
compensation amplification factor is e— (0.1 to 1.27). 

To keep £ = 0.707, when the ratio of the less significant 
to the most significant pole of the plant transfer function, 
r = pi/p2, varies from 50 to 1, the value of s may vary 
from 0.1 to 1.27. The case r = 2.5, s= 1 can be identified 
as shown in Figure 15 [23], [24], [25]. 

If the ratio is r = 2.5, as in the discussed case, then s- 1. 
If r < 2.5, ^is within the limits s- (1 to 1.27). If r >2.5, 
then s should be within the limits s- (0.1 to 1). 





Figure 15: Relationship between the Damping Ratio C, and the Factor s 
for Different Poles Ratios r = pi/p 2 

For the plant with the transfer function represented by 
Equation (26), a two- stage lead compensation and 
attenuation is applied. The poles pi = -50 and p2 = -20 
are cancelled and amplification factors ai = 0 L 2 = 10 and 
s=l are used. 

By applying Rule 1 and using a related sequence, the 
two-stage lead stage is presented by: 

' (1 + a,Ts)(\ + a^s) 

GJs) = - UJ. = 

a a (1 + 7[s)(l + T 2 s) 

(1 + 0.02^(1 + 0.05^) 

~ 100(l + 0.002^)(l + 0.005) 

Further from Rule 5: 

G (s) = saa =1x10x10 = 100 

c 1 2 


(27) 


(28) 


s(l + 0.02s)(l + 0.05s) 

Step 3: Optimization of £ , t s /t m and PMO 


For optimization of a marginal closed-loop system of 
Type 1 , Rule 1 is used in the same manner, while Rule 2 
is modified to Rule 5 and Rule 3 is ignored. 


Rule 1: To optimize t s /t m and the PMO of a Type 1 
marginal closed-loop system, a cascade multi-stage lead 
compensation with factors of 011,2, 3,... =10 should be 
applied for a zero-pole cancellation. The number of the 
compensating stages N should be one less than the order 
of the open-loop system N = n +m - 1 . 


The transfer function of the complete compensator as a 
product of the three sections is: 


G (s) = G ( s)xG (s) = 

c c c 

100(l + 0.02^)(l + 0.05^) 
“ 100(1 + 0.0025)(1 + 0.005s) 
(1 + 0.02^(1 + 0.05^) 

" (l + 0.002^)(l + 0.005^) 


(29) 


By applying the full compensation, considering Equations 
(26) and (29), the continuous prototype of the transfer 
function of the open-loop compensated system becomes: 
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GO) = G 0) x G 0 ) = 

c p 4 

70 


5(1 + 0.002s)(l + 0.005s) 


(30) 


Step 4: Conversion of the Continuous Prototype 
Compensated System into its Digital Equivalent. 

The transient responses of the closed-loop digital control 
system before and after the compensation are shown in 
Figure 16 and are obtained by the code: 

» Gp=50000/(s+50)/(s+20)/(s+0) 

» Gpfbd = c2d(Gpfb, 0.004, "tustin ) 

» G= 7000000/(s+500)/(s+200)/(s+0) 

» Gpfbd = c2d(Gpfb,0. 002, f tustin ) 

» step(Gpfbd,Gfbd) 


Step Responseof the Original and the 
Compensated Digital Control System 



Figure 16: Comparison between the Transient Responses of the 
Original and Compensated Digital System Type 1 

The Bode diagrams, after applying the compensation 
technique, are achieved by the code as shown below and 
are presented in Figure 17. 

» Gpd=c2d(Gp, 0.005, "\ tustin ) 

»Gd=c2d(G, 0.005, tustin ) 

» margin(Gd) 


Bode Diagram 


Gm = 20 dB (at 268 rad/sec) , Pm = 64.3 deg (at 65.3 rad/sec) 



Frequency (rad/sec) 


Figure 17: Bode Diagrams of the Digital Compensated System Proving 
Phase Margin PM = 64.3 and Damping £ « 0.643 


Comparison between the results and the objectives for 
optimal performance is summarized in Table 3. It is seen 
from Figure 16, Figure 17 and Table 3 that the results 
after the digital compensation are better or a close match 
with the objectives. 


Table 3 

Objectives and Real Results for Case 3 


Specifications 

Objectives 

Real Results 

Consideration 

c 

= 0.707 

= 0.643 

Close 

PMO 

<4% 

= 3.3% 

Better 

ts(l%)/tm 

< 1.49 

= 1.4 

Better 

e ss (t) 

= 0% 

= 0% 

Matching 


Step 5: Design of a Digital Compensation Controller 
based on a Microcontroller (Case 3). 

If considering Equation (29), the transfer function of the 
modified compensation controller stage can be presented 
as: 

_ (1 + 0.02^(1 + 0.05^) __ 

C ( S - (1 + 0.0025)(1 + 0.005s) “ 

O.OOls 2 +0.075 + 1 
0.0000 Is 2 + 0.007 + 1 

Further, the digital equivalent of the compensator 
controller stage is realized by the code: 


» Gc=tf([0.001 0.07 1], [0.00001 0.0071]) 
Transfer function: 

0.001 s A 2 + 0.07 s + 1 

le-005 s A 2 + 0.007 s + 1 
» Gcd = c2d(Gc, 0.001, ' tustin ) 

Transfer function: 

75.29 z A 2- 145.4 z + 70.2 

z A 2- 1.418 z + 0.4909 
Sampling time: 0.001 


Then, the transfer function of the series controller stage 
in the discrete-time domain is as follows: 

q 75.29z 2 - 145 ,4z + 70,2 _ 7(z) (32) 

z 2 -1.418z +0.4909 

Similarly to the other cases, a microcontroller can be 
connected in cascade with the original plant. To program 
the microcontroller, the transfer function of the series 
digital robust control stage is represented by the 
following difference equation [26], [27], [28]: 

y(kT) = -1.418y[(£ - l)r]+ 0.4909 y[(k - 2)7’] + 

+ 15. 29 x(kT) + ( 33 ) 

- 145 Ax[(k - 1)7’]+ 70.2x[(£ - 2)T ] 

Figure 18 demonstrates the analogue prototype of the 
compensated system. 



Figure 18: Block Diagram of the Analogue Prototype 
Compensated System 
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In Figure 19, the continuous model is converted into its 
digital equivalent by introducing the compensating 
microcontroller. 


This section describes the interactive procedure sequence 
[30], [31] for the analysis and design of the digital 
optimal compensation of the process control system. 


R(s) + 


A . 

Microcontroller 

Y(s) 

G P (s) 

C(s) 

y 




L 






Figure 19: Block Diagram of the Digital Compensated System 


4. Flowchart Diagram of the Analysis and 
Design Procedure Sequence 

The sequence of the analysis and the design procedure 
are described by the following flowchart diagram [29]. 



Figure 20: Flowchart Diagram of the Analysis and 
Compensator Design Procedure Sequence 


The system identification is an important element of any 
process of analysis. The dynamics of the plant process 
control system is described by a differential equation and 
further its transfer function Gp(s) can be derived. The 
ITAE criterion is considered for the design of the optimal 
compensator. 

Further, the system is identified as Type 0 or Type 1. For 
any one of the cases evaluation of the system is 
performed in terms of transient response, marginal gain 
and relative damping. The design procedure sequence 
applying Rule 1 and Rule 2 follows in terms of relative 
damping optimization and proper amplification. 

The constants a, y , S and s are determined from their 
graphical relationships with the DR, PMO and t s /t m . 
These constants are further used for the implementation 
of Rule 1, Rule 3, Rule 4 and Rule 5. 

The results are presented graphically after the initial 
system optimization, applying Step 1, Step 2 and Step 3 
and displaying the original system transient response and 
gain margin original with the aid of the D-partitioning 
stability analysis. 

If the system is Type 0, it is examined in terms of the 
difference between its dominant and its insignificant 
poles. If the difference is large, single stage lag 
compensation is applied. If the difference is insignificant, 
two-stage lag compensation is applied. For a system of 
Type 1, two-stage lead compensation is applied. 

Finally, considering the Euler's approximation for choice 
of the proper sampling period T s and applying the 
Bilinear Tustin Transform, digital conversion from the 
continuous-time presentation into the digital equivalent is 
performed by applying Step 4 and Step 5. 

Taking into account the compensator’s digital transfer 
function, its corresponding difference equation is derived 
based on which a microcontroller can be programmed to 
operate as a digital optimal compensator. 

The results after the digital conversion are presented 
graphically in the digital time-domain in terms of the 
system’s transient response and Bode diagram. 

Comparison is completed for the system’s performance 
before and after applying of the digital compensator and 
it is observed that an optimal compensation is 
accomplished. 

5. Conclusion 

The originality of the suggested technique of multi-stage 
compensation is based on the statement of a number of 
rules, which are applied in a predetermined sequence to 
some known theoretical procedures, like lead-lag 
compensation and zero-pole cancellation. 

By analyzing, combining and applying these rules to 
control systems of different nature, some new ideas are 
developed, bringing exceptional final results of the 
system’s performance. 
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The analogue compensation equipment consists of three 
major parts, facilitating the initial stage of the design. 

The analogue lead compensator section eliminates all 
insignificant poles of the original plant transfer function 
and introduces new properly designed dominant poles, 
improving the transient response and especially the 
relative damping ratio of the system. 

The amplifying section of the compensator reduces 
considerably the steady- state error. 

The analogue lag compensator section eliminates the 
most significant pole of the original plant transfer 
function and improves further the transient response. 

To upgrade further the compensator equipment, its 
analogue continuous transfer function is converted to its 
digital equivalent. This conversion is taking into account 
the Euler's approximation and is achieved by applying 
the Bilinear Tustin Transform, considered as one of the 
most accurate techniques. 

The excellent final results, observed from the digital 
system performance, prove once again the perfect match 
between the continuous compensator and its digital 
equivalent, due to the application of the Euler's 
technique. 

Difference equations are derived based on the digital 
equivalents of the compensators transfer functions. 
Microcontroller connected in cascade with the original 
plant can be programmed, based on the difference 
equation, reflecting the operation of the digital 
compensator. 

The latest manufactured microcontrollers offer very high 
operating frequencies that can handle successfully the 
operation of each one of the discussed in this research 
control systems. For example the TMS320F2837xS is a 
powerful 32-bit floating-point microcontroller designed 
for advanced closed-loop control applications that 
provides 200 MHz of signal processing performance [32] 
and can be used for each one of the considered cases. 

There is significant number of advantages associated 
with the digital control systems, compared with the 
continuous control systems. Due to this reason, the 
choice of compensation design strategy, resulting in the 
implementation of microcontrollers, is the better and 
preferable option and can be applied for any process 
control system. 
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